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I/"") ■ A non-Hermitean operator does not necessarily have a complete set of 

eigenstates, contrary to a Hermitean one. An algorithm is presented which 
| allows one to decide whether the eigenstates of a given PT-invariant operator 

on a finite-dimensional space are complete or not. In other words, the 
algorithm checks whether a given PT-symmetric matrix is diagonalizable. 
The procedure neither requires to calculate any single eigenvalue nor any 
numerical approximation. 
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I. Introduction 
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The physical interpretation of PT-invariant operators — and hence their relevance for the 
description of physical systems — continues to be debated [II 121 El- There is, however, no 
doubt about the cathartic role of PT-symmetry: it has become more evident what it means 
to let go hermiticity in exchange for a weaker property such as PT-invariance. The suc- 
cess and ease to describe quantum mechanical systems in terms of hermitean operators is 
based on two of their generic properties, namely the existence of real eigenvalues and their 



diagonalizability, i.e. the completeness of their orthonormal eigenstates. These properties 
do not persist if a quantum system was described by a PT-symmetric Hamiltonian: its 
eigenvalues could be complex, and its eigenfunctions would, in general, neither be pair- 
^ ■ wise orthogonal nor form a complete set. Given a PT-invariant operator, it thus appears 

desirable to decide whether it is diagonalizable or not. 

The purpose of this contribution is to provide an algorithm answering the question 
of whether a given PT-invariant Hamiltonian operator in a finite-dimensional space does 
or does not possess a complete set of eigenstates. It is convenient to represent such an 
operator as a PT-symmetric matrix M, say. A procedure will be outlined which, after a 
finite number of steps, will announce whether the matrix M at hand is diagonalizable or 
not. In principle, the algorithm can be carried out by hand for matrices of any dimension, 
and no approximations are necessary. 

Often, the question of diagonalizability will arise in a more general setting where one 
considers not just a single matrix but a family of PT-symmetric matrices M(e),e £ JR.. The 
parameter e measures the strength e £ K. of a "perturbation" which destroys hermiticty 
while repecting PT-invariance. As the parameter varies, all of the cases described previ- 
ously may occur: typically, two real eigenvalues merge into a single real one at a critical 
value of e, subsequently splitting into a pair of two complex conjugate eigenvalues, or vice 
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versa. These dramatic modifications are accompanied by changes in the nature of the 
eigenstates of the PT-invariant operator, possibly no longer spanning the space on which 
M(e) acts. 

This behaviour can be understood in terms of so-called exceptional points j3] which are 
known to occur when a matrix is subjected to the perturbation depending analytically on 
a parameter such as e. At such a point, the corresponding matrix is not diagonalizable, 
and its spectrum may undergo a qualitative change. For a hermitean operator subjected 
to a parameter-dependent hermitean perturbation, exceptional points cannot occur. 

If one applies the algorithm testing for diagonalizability to a parameter-dependent ma- 
trix M(e), it will output a polynomial in e instead of a number. Its zeros correspond to 
those values of the perturbation parameter where the matrix family M (e) has exceptional 
points. The matrices corresponding to these values of the perturbation are not diagonal- 
izable, and the spectra of matrices for nearby values of the parameter differ qualitatively. 

The following section summarizes the properties of PT-invariant systems in terms of 
(2 x 2) matrices. Then, the link between diagonalizability and the so-called minimal 
polynomial is reviewed. In Section 3, the algorithmic test is presented which consists of 
constructing the minimal polynomial of the matrix followed by a search for degenerate 
roots by means of the Euclidean algorithm. Various methods are known to effectively 
calculate the minimal polynomial of a matrix, outlined in Section 4. Simple examples 
are studied in Section 5, leading to some general conclusions about the structure of PT- 
symmetric Hamiltonian operators in finite-dimensional spaces. Section 6 summarizes the 
results and discusses the challenge to extend them to state spaces of infinite dimension. 

II. PT-invariant systems 

A matrix H is PT-invariant 

[H,PT] = 0, (1) 

if it commutes with the product of parity P and the anti-unitary operation of time re- 
versal T, represented here by complex conjugation, T^TH = H*. Eq. implies that 
the characteristic polynomial of any PT-symmetric operator H has real coefficients only. 
Consequently, its roots are either real or come in complex-conjugate pairs. One way to 
show this is to construct a basis in which the Hamiltonian has real matrix elements only 

Let us briefly review the properties of PT-symmetric systems by considering the most 
general PT-invariant matrix of dimension 2, 

with parity given by the Pauli matrix a x in the standard representation. For real numbers 
a and b, the matrix H is not only PT-invariant but also hermitean. Thus, its eigenvalues 
are real, and its orthonormal eigenstates span C 2 . For a* ^ a and b = 0, H has a pair of 
complex conjugate eigenvalues and two orthonormal eigenstates. Matrices of the form 

H=( j _ & .) , be[-l,l], (3) 
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are particularly interesting. For |6| < 1, one finds a pair of two complex conjugate eigen- 
values, 

E± = ±Vb 2 - 1 e z'R , (4) 
associated with two non-orthogonal eigenstates, 



b \ 1/6 
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(5) 



(6) 



When 6 = ±1 in Q, H has a two-fold degenerate eigenvalue, E = 0, and there is only 
one eigenstate, namely, 

-( Tl 

This situation, impossible for a hermitean matrix, is usually described by saying that the 
algebraic multiplicity of the eigenvalue Eq is two while its geometric multiplicity equals one: 
the characteristic polynomial of M has a double root associated with a single eigenvector 
only. In this case, the matrix H is not diagonalizable: a similarity transformation sending 
it to a diagonal matrix cannot exist since its eigenstates would span the space C 2 . 

III. Diagonaliza bility and the minimal polynomial of a matrix 

Each square matrix M of dimension N satisfies the identity 

Pm(M) = 0, (7) 

where £>m(A) is the characteristic polynomial of M, 

Pm(A) = det (AE — M) , (8) 

with E being the unit matrix of dimension N. In other words, the characteristic poly- 
nomial of M annihilates the matrix M. The polynomial Pm(A) has degree N and it is a 
monic polynomial, that is, the coefficient multiplying the highest power of A is equal to 1. 
Obviously, many other monic polynomials of higher degree also annihilate M: simply take 
Pivi(A),pivi(A), • • • K is less obvious, however, whether one can find polynomials of degree 
less than N which annihilate M. This, in fact, depends on the properties of the matrix M. 

Define [7] the minimal polynomial of the matrix M as the monic polynomial ?7Jm(A) of 
least degree which annihilates M : 

m M (M)=0. (9) 

The minimal polynomial my\(X) is unique [JJ, and its degree N is less than or equal to 
the degree of the characteristic polynomial, iVo < N. The minimal polynomial divides the 
characteristic polynomial without a remainder, m M (A) |pm(A), or equivalently, 

Pm(X) = rf M (A)m M (A) , (10) 

where c?m(A) is a non-zero polynomial of degree less than N. The characteristic and the 
minimal polynomial of the matrix M coincide in the case where g?m(A) = 1. 
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In general, the minimal polynomial has u roots M u , 

m M {\) = \[{\-M v Y% u <N, (11) 

with multiplicities fi u summing to N = /ii +/i 2 + • • • + AW Here is the important property 
of the polynomial 771m(A): the matrix M is diagonalizable if and only if each root M v in 
fITTj) has multiplicity one, //„ = 1, v — 1, . . . , i/ , that is, 

m M (A) = JJ(A - M„) , all M v distinct . (12) 

«/=i 

No polynomial of degree less than tom(A) annihilates the matrix M. 

Let us illustrate the properties of minimal polynomials using low-dimensional matrices. 
Consider the matrix A with entries (1,1,2) on the diagonal, and zero elsewhere. Its 
characteristic polynomial is given by 

p A (A) = (A-l) 2 (A-2), (13) 

while its minimal polynomial reads 

m A (A) = (A-l)(A-2), (14) 

being of the form (fTTJl . with Nq = v$ = 2. This is easy to verify since m A (A) = A 2 — 3 A + 
2E = holds, while none of its factors annihilates M: both (M — E) and (M — 2E) are 
different from zero. Thus, the minimal polynomial divides the characteristic one, p A (A) = 
(A — l)m A (A), leading to d A (A) = (A — 1). Due to (JHJ), the matrix A is diagonalizable — a 
correct but hardly surprising result since the matrix A has been diagonal from the outset. 
Here is the instructive part of the example: consider the matrix 



B= 1 , beC, (15) 




which is different from A as long as b is different from zero. The characteristic polynomial 
of B equals that of A but the matrix B must have a different minimal polynomial since 
m A (B) 7^ 0. It is not difficult to verify that no linear or quadratic polynomial annihilates B 
as long as b ^ 0. This implies that its minimal polynomial coincides with its characteristic 
polynomial, 

p B (A) = m B (A), cf B (A) = l. (16) 

Consequently, the minimal polynomial of B does not have the form specified in (jl2j) . and 
the matrix B is not similar to a diagonal matrix. Inspection shows that B indeed contains 
a (2 x 2) Jordan block for any nonzero value of b. 

For PT-invariant matrices, both the polynomials o?h(A) and m^(X) have real coefficients 
only, just as the characteristic polynomial. This will be shown once the function g?m(A) in 
(fTU|) has been defined in general (cf. Sec. IV. 2). 
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IV. An algorithmic test for diagonalizability 



A square matrix M of dimension N is diagonalizable if its minimal polynomial is a product 
of factors (A — M v ) with all numbers M u , v — 1, 2, . . . , n < Vq, distinct, as shown in Eq. 
()12|) . Consequently, to test for diagonalizability of given a matrix M, one needs to 

(2) find the minimal polynomial mivi(A) of the matrix M; 

{11) determine whether the polynomial mivi(A) has single roots only. 

To calculate numerically the roots of either the characteristic or the minimal polynomial is 
not a valid approach since, in general, the exact roots of a polynomial cannot be specified 
in a finite procedure. Any algorithmic implementation must generate answers to (1) and 
{ii) in a finite number of steps. Note that even if the first step has been implemented, it 
is unlikely that the minimal polynomial will emerge in factorized form. 

Interestingly, it is possible to construct the minimal polynomial of a matrix and to 
check for degenerate roots in a finite number of steps. In both cases one searches for 
common factors of polynomials, which is achieved algorithmically by the Euclidean division 
algorithm. These results seem to have been put together for the first time in jS] in order to 
decide algorithmically whether a given matrix is diagonalizable. As it stands, it could be 
applied to the non-hermitean matrix govering the motion of two coupled damped classical 
oscillators studied in [§]. 

In the following, a slightly simplified approach to the problem of diagonalizability is 
presented, adapted to matrices with PT-symmetry. Before implementing the steps {%) 
and (tt), the Euclidean algorithm for polynomials will be presented briefly to establish 
notation. 

IV. 1 Euclidean division algorithm for polynomials 

Given two integer numbers po > pi, say, the Euclidean division algorithm outputs their 
greatest common divisor, denoted by gcd(po,Pi) £ No, after a finite number of steps. It 
works as follows: first, you need to express the larger number as gi-fold multiple of the 
smaller number plus a remainder p 2 , 

Po = qiPi+P2, <?i,P2eN , Pi>p 2 <0- (17) 

This relation implies that any common divisor of po and p\ divides p 2 as well, hence 
gcd(p ,Pi) = gcd(pi,p 2 ). Thus, it is sufficient to search for the greatest common divisor 
of the pair {pi,P2)- This can be achieved by increasing each index in (|17|) and feeding in 
the pair (^1,^2) instead of (po,Pi), etc. Since p > p\ and p 1 > p 2 , the algorithm will 
stop after a finite number of iterations and produce a remainder equal to zero, pt+i = 0, 
say. The non-zero remainder pk generated in the penultimate step is the desired result, 
gcd(po,Pi) = Pk- If gcd(p ,Pi) = 1, the numbers p and p\ are relatively prime, otherwise 
a common divisor different from one has been identified. 

A polynomial in the variable A can be written as a unique product of linear factors (A — 
A„) where the numbers A„ G C are its roots. This representation makes polynomials similar 
to integer numbers in some respects. The equivalent of the Euclidean algorithm, when 
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applied to two polynomials, outputs their greatest common divisor, which is a polynomial 
itself. This result is based on the fact that any two polynomials po(\) and pi(A), with 
degpo(A) > degpi(A) are related by 

Po(A) = gi(A)pi(A) +p 2 (A) , deg Pl (A) > degp 2 (A) > , (18) 

which is the equivalent of (fTTj) . The polynomials qi(X), with deggi(A) = (degpo(A) — 
degpi(A)), and hence p 2 (A), are found from long division. If po(X) and pi(A) have a 
common factor, then p 2 (A) must have this factor as well. Thus, it is sufficient to search for 
gcd(pi(A),p2(A)) instead of gcd(p (A),pi(A)) but the degrees of the polynomials involved 
have effectively been reduced. Consequently, this procedure can be repeated all over again 
and it halts once a vanishing remainder has been obtained, pk+i{X) = 0, say. Then, the 
greatest common factor of the polynomials po(\) and pi(A) is given by the last non-zero 
remainder polynomial, Pk{X), calculated in the next-to-last application of the algorithm. 
If degpfc(A) = 0, the initial polynomials are "relatively prime," otherwise their greatest 
common divisor is a polynomial of degree at least one. 



IV. 2 Step (z): Finding the minimal polynomial of a matrix 

The function g?m(A) relates the minimal polynomial wim(A) of the matrix M to its character- 
istic polynomial pu{X) according to Eq. (jlOj) . Hence, the minimal polynomial associated 
with M is known once the characteristic polynomial and the function g?m(A) have been 
determined. 

Two steps are required to construct the function g?m(A) [7j. First, you need to calculate 
the matrix Dm = adj (AE — M), given by the transposed cofactors — or signed minors — of 
the matrix (AE — M). The adjoint of a matrix, C say, always exists, and it satisfies the 
relation 

Cadj C = (detC)E. (19) 

For det C 7^ 0, Eq. (|19|) leads to the familiar expression of the inverse matrix of C. 

According to [7] the polynomial du{X) is given by the greatest (monic) common divisor 
of the N 2 elements of adj (AE — M), 

d M (X) =gcd{(D M )nmhm = l,...,N} , (20) 

Thus, in a second step, you need to apply the Euclidean algorithm to all pairs of entries 
of the matrix adj (AE — M). Having thus identified the function g?m(A), the minimum 
polynomial of M follows from (jl(Jj) . 

mM{x) = MX)- (21) 

Now it is possible to show that, for a PT-invariant matrix, the polynomials c£m(A) and 
?7Jm(A) have real coefficients only, just as the characteristic polynomial. Using a basis in 
which all elements of H are real, leads to 

(D H (A))* = (adj(AE - H))* = adj(A*E - H) = D H (A*) . (22) 
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which states that in this basis the adjoint of H has only real matrix elements (except for 
the unknown A). Taking (|20|) into account this leads to 



(d H (A))* = gcd {(D H ) nm (A*)|n, m = l,...,N} = d H (X*) , (23) 

which, in conjunction with (p H (X))* = Ph(\*) and Eq. (|2*T|) implies indeed (m H (A))* = 
m H (A*). 

Let us verify that this procedure outputs the correct minimal polynomials for the 
matrices A and B introduced in Section 3. The adjoint of the matrix (AE — B) reads 

/ A - 1 -b 
adj (AE - B) = adj A - 1 

\ A-2 

(A-l)(A-2) 6(A-2) 

(A-l)(A-2) j . (24) 

(A-l) 2 

Due to the simplicity of the matrices involved, the Euclidean algorithm can be run "by 
inspection:" for 6^0, the only common factor among the entries in (}2"4*|) is given by 
^b(A) = 1. Consequently, the minimal and the characteristic polynomial of B coincide 
as stated in Eq. (jlfij) . If the parameter b takes the value zero, B turns into A, and and 
a non-constant greatest common divisor emerges, g?a(A) = (A — 1). Using Eq. (JTUJ), one 
obtains the minimal polynomial m^(X) = (A — 1)(A — 2), agreeing with Eq. (|14|). 

In |8j, a different approach to determine the minimal polynomial of a matrix M has 
been presented which, ultimately, is also based on finding the greatest common divisor of 
specific polynomials. According to ^U], any method to determine whether the matrices 
M° = E, M, M 2 , . . . , M^ -1 , are linearly dependent, can be used to construct the minimal 
polynomial of M; two such methods are described in this reference, and a third one can 
be found in The latter approaches have in common that they are not based on the 
Euclidean algorithm. For actual calculations, it is convenient to resort to a Mathematica 
program 12] to find the minimal polynomial of a matrix M. 

4.3 Step (m): Identifying degenerate roots of a polynomial 

Once the minimal polynomial 771m(A) has been found, one needs an algorithm to decide 
whether it has single roots only jH]. Imagine a polynomial m(A) to have an s-fold root A , 
2 < s < N. Its factorization reads 

m(A) = (A-A ) s ... , (25) 

where the dots indicate a polynomial of degree (N — s). Its derivative takes the form 

^ = (A-A r 1 ..., (26) 

the dots standing again for some polynomial of degree (N — s) . Obviously, the polynomial 
and its derivative are not relatively prime: m(A) and ra'(A) have a factor (A — Aq) 5 " 1 
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of order at least one in common. Thus, applying the division algorithm to the pair 
(mivi(A), m' M (X)) checks whether the polynomial mivi(A) has the form (fT2j) . If the pro- 
cedure outputs gcd(mM(A), m' M (A)) oc 1, all roots of are distinct and the associated 
matrix M is diagonalizable, otherwise it is not. 

This concludes the description of an algorithm to test for diagonalizability of a given 
PT-symmetric matrix M. No fundamental changes are necessary if one studies a parameter- 
dependent family of matrices M(e). However, the algorithm will output conditions poly- 
nomial in the parameter e, indicating specific parameter values where diagonalizability 
breaks down. It is convenient to study the resulting modifications by working out some 
simple examples, illustrating at the same time the proposed algorithm. 



V. Examples 

V.l Matrices of dimension (2 x 2) 

Let us apply the algorithm described above to the matrix H in (0) assuming the numbers 
a and b to be different from zero. Its characteristic polynomial reads 

j9 H (A) = det(AE - H) = A 2 - 2(3?a)A + |a| 2 - \b\ 2 , (27) 

while its minimal polynomial is found via the function dy\(X) equal to the highest common 
factor of the matrix 

D„(A)=ad j( AE-H,= ( A _-;* \~— a 

By inspection, a non-constant factor only exists among the four entries of Dh(A) if b = 
Qa = 0. In this case, H turns into a real multiple of the identity, hence it is diagonalizable. 
This observation illustrates a fine point of the construction of the minimal polynomial: 
even upon identifying a non-constant function g?h(A), the minimal polynomial ?tih(A) may 
still be of the form (|12|) . For b = = 0, the characteristic polynomial turns into (A— 9?a) 2 , 
implying indeed m^(X) = (A — 3?a). Here, the function d^(X) = (A — 3?a) removes factors 
of the characteristic polynomial which stem from the degeneracy of an eigenvalue of H. 

From now on, either b or are assumed to be different from zero, hence c?h(A) = 1, 
and the minimal polynomial my\(X) is given by Eq. (j2*7|) . 

m H (A)=p H (A), ( 29 ) 

which concludes the first step of the algorithm. 

In the second step of the algorithm, the search for multiple roots of mn(A) = po(X), 
one needs to determine the highest common factor of the minimal polynomial and its 
derivative, 

m' H (A) = 2(A — 3?a) = Pi(X) . (30) 

Applying the Euclidean algorithm to the pair (m^(X) , m' H (X)) means to solve for a poly- 
nomial qi(X) = AX + B and for p2(A) such that 

m H (A) = (AX + B)m' H (X)+p 2 (X), A,BeR. (31) 
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The unknowns are easily obtained as 



A=i, B = ~sRa, p 2 (A) = (3a) 2 - |6| 2 ~ A . (32) 

Two possibilities now arise: either p 2 (A) equals zero or it does not. The first case occurs 
if 

(3a) 2 = |6| 2 , (33) 

and the algorithm comes to a halt. As mentioned above, the greatest common divisor of 
the initial polynomials is then given by the penultimate (monic) remainder polynomial, 
i.e. pi(A) = (A — dta). It follows that toh(A) and its derivative do have a common non- 
constant divisor, so that H cannot be brought to diagonal form. It is easy to verify that 
mn(A) = (A — 3?a) 2 when (}3*3*j) holds, confirming that the minimal polynomial of H has a 
double root dta. Furthermore, a simple calculation shows that the matrix H has indeed 
a single eigenstate only if the relation = ±|6| holds. Note that this result covers 
the example of a non-diagonalizable H of Section 2, where a = i and b = ±1 had been 
considered. 

Finally, if a and b do not satisfy ([33)1. the remainder polynomial p 2 (A) does not vanish. 
Being a constant, the algorithm is bound to stop after the next iteration. Determine 
02(A) = (CA + D) and pz{\) such that 

(A-SRa) = (CA + L>)((3a) 2 - |6| 2 ) +p 3 {\), C,DeR, (34) 

holds, i.e., 

1 — 9£a 

C = (3a)* - |&r D= (3a)'-|6|" P3(A) = °' (35) 

The algorithm halts indeed due to P3(A) = 0, and the penultimate remainder is p 2 (A) oc 1, 
indicating that the minimal polynomial does not have any degenerate roots, and H is 
diagonalizable. 

In summary, the matrix H is diagonalizable for all parameter values except when 3a = 
±|6|. In this case the algebraic multiplicity of its eigenvalue is two, while its geometric 
multiplicity is one; otherwise the multiplicities both equal two. It is important to note 
that it was not necessary at any stage to determine the eigenvalues of H. 



V.2 Matrices of dimension (4 x 4) 

It is instructive to apply the algorithm to a PT-invariant matrix of dimension 4, 



H = 
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which depends on a perturbation parameter e. As before, the action of T on a matrix 
effects complex conjugation of its entries, while P is now given by a (4 x 4) matrix with 
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entries equal to one along its minor diagonal and zero elsewhere. Eq. (JT} is then readily 
verified. The characteristic polynomial of H is given by 



p H (A) = A 4 + a £ X 2 + (3 £ , (37) 

with a £ = 2e 2 — 2s 2 — 5 2 and f3 E = e 4 — (2s 2 + 5 2 )e 2 + s 4 . The minimal polynomial of H 
coincides with the characteristic one, toh(A) = Ph{X), since the only common factor of the 
matrix elements of Dh is equal to one, eta (A) = 1. To see this, it is sufficient to calculate 
the two matrix elements [D H ]i4 = —s 2 5, and [D H ]23 = — (X 2 + e 2 )S, for example. Whatever 
the value of e, for nonzero s and 8 the only common divisor is one, so that eta (A) = 1. 

Let us now determine gcd(mi-i(A), m' H (X)) by the Euclidean algorithm, where my\(X) = 
Po(A) is given in Eq. ([3*7)1 and m' H (X) = Pi(X) = 4A 3 + 2a £ X. Comparing powers of A in 
Eq. ("21} with the polynomials just defined, one obtains 

A = \, B = 0, P2 (X) = ^-X 2 + (3 £ . (38) 

The algorithm only stops here if a £ = f3 £ = which would require s = 5 = 0, contrary to 
both s and 5 being different from zero. If a £ = is assumed, H is diagonalizable for all e 
since f3 e cannot take the value zero, and the algorithm stops after the next step, outputting 
f3 £ oc 1 as greatest common factor. Assume now a £ ^ and apply the division algorithm 
to the pair (pi{X) , p 2 (X)) . The unknown constants in (72(A) = CA + D, and the remainder 
polynomial P3(A) are found to be 

C=—, D = 0, p 3 (X) = —(a 2 £ -4P £ )X. (39) 

a £ a £ 

For the algorithm to stop, one must have ps(X) = 0. This, however, does not happen 
whatever the value of e since a 2 — 4j3 £ = 5 2 (As 2 + 5 2 ) > 0. The next iteration of the 
algorithm leads to 

E= wrw Y F = ' ^ (A) = ft - (40) 

where 93(A) = EX + F. Producing a remainder polynomial of degree zero in A, the 
condition for the minimal polynomial to have multiple roots is finally given by 

P e = e 4 - (2s 2 + 5 2 )e 2 + s 4 = . (41) 

This fourth-order polynomial in e has roots 

±e± = ±Vv± Act, a = s 2 + — >0, Aa = Va 2 - s 4 e (0, a) . (42) 

For each of these four real values of the parameter e, the matrix H is not diagonalizable. 
In other words, the algebraic and geometric multiplicity of the eigenvalues of H do not 
coincide, and its eigenstates do not span the space C 4 . It is important to note that the 
eigenvalues of the matrix H are not known at this stage. 
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Figure 1: Parameter Space of the matrix H denned in (|3fi|) : region I: four real eigenvalues; 
regions II: two real and one pair of complex-conjugate eigenvalues; region III: a two pairs 
of complex-conjugate eigenvalues; the matrix H is not diagonalizable on the full lines 
separating the regions I, II, and III 

V.3 The global structure of H 

Let us now determine the global properties of H. This is easily done upon combining (}4"2*|) 
with the characteristic polynomial ()37|) in its factorized form, 

p H (A) = (A - A + )(A + A + )(A - A_)(A + A_) (43) 

with roots 

±A± = ±y/a ±Aa-e 2 } (44) 

expressed directly in terms of a and Act. 

To graphically represent the parameter space of H and its properties, it is convenient 
to eliminate the parameter s by the scaling e — > se and 5 —>■ s5. This effectively amounts 
to sending s — > 1, and Eq. (jUJ) simplifies to e 4 — 2(1 + 5 2 /2)e 2 + 1 = 0, plotted in Fig. 1. 

Imagine to move along the dashed vertical line, determined by a fixed positive value 
of 5 > and variable e. For e = 0, the matrix H is hermitean, hence it has four distinct 
real eigenvalues and four orthonormal eigenstates. In region I, where < |e| < \[a — Act, 
Eq. (|44jl says that the eigenvalues remain real and distinct; a complete, not necessarily 
orthonormal set of four eigenstates continues to exist since (}4~T|) does not hold. When e = 
±y/ (a — Act), two eigenvalues coincide numerically, and the corresponding two eigenstates 
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merge into a single one, leaving H with an incomplete basis. Then, for y/a — Act < |e| < 
\/ a + Act, in region II, two real eigenvalues and a pair of complex-conjugate ones exist, 
with H being diagonalizable throughout since (|%TJl does not hold. At e — ±a/(ct + Act), 
the remaining two real eigenvalues degenerate to a single one, leaving H non-diagonalizable 
again with only three eigenstates. Finally, in region III, defined by a/ct + Act < |e|, the 
matrix H is diagonalizable and it comes with two pairs of complex-conjugate eigenvalues. 

Finally, for 5 = 0, the matrix H in (JHUj) decouples into an pair of identical two- 
dimensional matrices. The left boundary of region I sees the real eigenvalues of H de- 
generate pairwise which is consistent with the observations made earlier. At e = ±1, only 
two eigenstates exist while all four eigenvalues coincide numerically. Beyond this value of 
e, there are two pairs of identical complex-conjugate eigenvalues, and the associated basis 
is complete. 

For PT-symmetric systems described by matrices of higher dimensions it is, in general, 
not possible to find the roots of the characteristic polynomial. Nevertheless, a discussion 
of the parameter space can still be given: to this end one needs to detect the number of 
real and complex eigenvalues for each set of parameter values; an algorithm capable of 
doing this will be presented in 

VI. Discussion and Outlook 

An algorithm has been presented which allows one to determine whether a given PT- 
invariant matrix M is diagonalizable. To do so, it is not necessary to determine the roots 
of its characteristic polynomial. In terms of Linear Algebra, the algorithm decides whether 
the given matrix is similar [7j to a diagonal matrix or to a matrix containing at least one 
Jordan block of dimension two or more. Somewhat surprisingly, this question seems to 
have been addressed only recently from an algorithmic point of view. 

It seems worthwhile to point out that the test for multiple roots of a polynomials can, in 
fact, be used without any change to determine whether the eigenvalues of a given hermitean 
matrix are degenerate or not. The present author is not aware that this observation has 
been made before. 

When applied to a family of PT-symmetric matrices, the algorithm outputs polynomial 
conditions on the perturbation parameter. These conditions are satisfied for sets of matri- 
ces all of which are not diagonalizable, and they divide the full parameter space into regions 
of diagonalizable matrices with qualitatively different spectra. When combined with an 
algorithm to identify the number of real and complex eigenvalues of M, a complete picture 
of the system's properties in the entire parameter space can be established. 

Many PT-symmetric systems-including the first one studied from this perspective j^j- 
have been defined on Hilbert spaces with countably infinite dimension. Various concepts 
such as eigenvalues and eigenstates, or the difference between algebraic and geometric 
multiplicities of degenerate eigenvalues continue to exist in the more general case 0]. In 
spite of a close similarity of hermitean operators in finite- and infinite-dimensional spaces, 
many concepts of the matrix case are not easily carried over to the more general situation. 
For any algorithm, finiteness is a crucial feature: the number of steps required to identify a 
potential common factor of two polynomials is always finite, no matter what their degree. 
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It will be interesting to see whether algorithmic tests for diagonalizability of operators 
acting on spaces with countably infinite dimension can be found. 
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